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Lagrangian chaos and turbulence in fluid dynamics 


A. Bershadskii 
ICAR, P.O. Box 31155, Jerusalem 91000, Israel 


Randomization of the Lagrangian chaos in fluid dynamics has been analyzed using results of 
direct numerical simulations, laboratory measurements, and oceanic observations. The notion of 
distributed chaos has been used in order to quantify this phenomenon (the main dimensionless pa- 
rameter of distributed chaos has been used as a measure of randomization). The analysis includes 
isotropic homogeneous fluid motion, buoyancy-driven flows, particle-laden flows, magnetohydro- 
dynamic flows, and subsurface oceanic mixing. The role of the dissipative Birkhoff-Saffman and 
Loitsyanskii invariants as well as the invariants related to the spontaneous breaking of the local re- 
flectional symmetry (in particular, the Levich-Tsinober invariant) in the randomization phenomenon 


has been also investigated. 


I. INTRODUCTION 


Lagrangian approach can be very useful in studying 
mixing and dispersion in fluid dynamics. In a Lagrangian 
approach, an observer follows the motion (the trajecto- 
ries time evolution) of fluid particles, unlike an Eulerian 
approach where an observer considers the velocity field 
u(x,t) in a space domain at any time t. It is known 
that even a regular velocity field can produce rather 
irregular motion of the corresponding fluid particles (see, 
for instance, paper Ref. [a}). Sensitivity to the initial 
conditions (the exponential divergence of the initially 
nearby trajectories of the fluid particles) is the main 
characteristic of the Lagrangian chaos. 


The notion of smoothness can be used as a basis for 
classification of the chaotic regimes. Smooth dynamics 
is generally characterized by the stretched exponential 
Lagrangian frequency spectra 


E(w) x exp —(w/we)?. (1) 


where 1 > 6 > O and w is the frequency. The value 
8 =1 (i.e. the exponential spectrum): 


E(w) x exp(—w/w-), (2) 


can be associated with deterministic chaos (see, for in- 
stance, Refs. (2]--[5}). 

For 8 < 1 the Lagrangian chaotic-like dynamics is not 
deterministic but nevertheless smooth (the distributed 
chaos, see below). The non-smooth randomized dynam- 
ics is usually characterized by the scaling (power-law) 
spectra (hard turbulence [6]). 


In the distributed chaos (ie. for G < 1) the value 
of @ can characterize a measure of randomization. 
Namely, lesser values of the dimensionless parameter 
8 correspond to stronger randomization of the fluid 
dynamics. One can see in Fig. 1 that smaller values 
of the Taylor-Reynolds number Re, correspond to 
larger values of 8 up to the 6 = 1 (corresponding to 
deterministic chaos). It seems rather natural that larger 
values of R, should generally correspond to stronger 
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FIG. 1: Lagrangian frequency spectra of velocity in statisti- 
cally stationary isotropic homogeneous fluid motion (DNS at 
different values of Ry). The spectra are vertically shifted for 
clarity. 


randomization. 


Figure 1 shows the results of direct numerical sim- 
ulations (DNS) of statistically stationary isotropic 
homogeneous chaotic/turbulent motion of an_ incom- 
pressible fluid reported in the paper Ref. . The 
spectral data were taken from Fig. 5 of the Ref. (T, is 
the Kolmogorov time scale). The DNS was performed for 
the Navier-Stokes equations. The statistically stationary 
chaotic/turbulent motion was maintained in the range 
38 < Re, < 234 using a large-scale forcing. The dashed 
curves in Fig. 1 indicate the stretched exponential 
spectra Eq. (1) and the exponential spectrum Eq. 
(2) (6 = 1, the deterministic chaos for the smallest 
considered Re, = 38). 


In the next sections, the values of 8 shown in the Fig. 1 
will be related to different invariants of the fluid dynamics 
and it will be also shown that spontaneous breaking of 
the local reflectional symmetry plays an important role in 
the fluid dynamics (including the isotropic homogeneous 
motions for sufficiently large Re). 
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FIG. 2: Lagrangian frequency spectrum of zonal and merid- 
ional velocity for trajectories obtained using the subsurface 
floats at a depth between 400 and 1000 m in the Atlantic 
Ocean (the World Ocean Experiment - WOCE). 


II. LAGRANGIAN DETERMINISTIC CHAOS 
IN MORE COMPLEX FLOWS 


It was shown in the Introduction that the Lagrangian 
deterministic chaos (the exponential spectrum) takes 
place in isotropic homogeneous fluid motion at small Re 
(about the Eulerian wavenumber exponential spectra in 
the isotropic homogeneous fluid motion at small Re 
see, for instance, Ref. [8}). 


Let us now give several other examples of the exponen- 
tial Lagrangian frequency spectra in more complex flows. 


Figure 2 shows the Lagrangian frequency spectrum of 
zonal and meridional velocity for trajectories obtained 
using the subsurface floats at a depth between 400 
and 1000 m in the Atlantic Ocean (the World Ocean 
Experiment - WOCE). The spectral data were taken 
from Fig. 3 of Ref. The floats were neutrally 
buoyant and ballasted for a specific pressure. A total 
of 132 736 daily observations were used to compute the 
spectra. 


The dashed curve in Fig. 2 indicates the exponential 
spectrum Eq. (2) (the deterministic chaos). The dotted 
arrow indicates the characteristic frequency we. 


Figure 3 shows the Lagrangian frequency spectra ob- 
tained in a direct numerical simulation (DNS) reported 
in a paper Ref. (Wy = 1/T). The spectral data 
were taken from Fig. 13b of the Ref. (10). 


In this DNS a particle-laden decaying isotropic 
homogeneous flow in a cubic box (with periodic bound- 
ary conditions) was studied using the incompressible 
Navier-Stokes equations (one-way coupling case) and the 
complete equation of particle motion (including pressure 
and viscous drag forces, inertia force of added mass, 
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FIG. 3: The Lagrangian frequency spectra computed for the 
velocity of fluid surrounding the inertial particles of corn and 
those of copper for a particle-laden decaying isotropic homo- 
geneous flow in a cubic box. The spectra are vertically shifted 
for clarity. 


force due to viscous stresses and fluid pressure gradient, 
viscous force related to unsteady relative acceleration, 
and gravity (or buoyancy) force). The initial Re, was 
equal to 25 and final Re, = 16. 


The Lagrangian frequency spectra shown in the Fig. 3 
were computed for the velocity of fluid surrounding the 
inertial particles of corn and those of copper (i.e. the 
local fluid existing around the solid point-like particle 
at each position of the particle along its trajectory). 
The dashed curves in Fig. 3 indicate the exponential 
spectrum Eq. (2) (the deterministic chaos). We will 
return to this DNS in Section IV. 


III. LAGRANGIAN DISTRIBUTED CHAOS 
AND DYNAMIC INVARIANTS 


When the characteristic frequency w. in the Eq. (2) 
randomly fluctuates one needs an ensemble averaging 


E(w) « [ P(w.) exp —(w/we-) dw (3) 


in order to obtain the average spectrum. 

Providing that the randomized dynamics is still 
smooth, the spectrum E(w) has the stretched exponen- 
tial form Eq. (1). Comparison of the Eq. (1) and Eq. 
(3) gives asymptote of the probability distribution P(w,) 
for large characteristic frequencies w, 


Plc) 0 hy +9/20-8)1 exp(—ykB/O-%), (4) 
where ¥ is a constant. 


The dissipative dynamics described by the Navier- 
Stokes equations has two fundamental invariants: 
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FIG. 4: The Lagrangian frequency spectra computed for the 
velocity of fluid surrounding the inertial particles of corn and 
those of uranium for a particle-laden decaying isotropic ho- 
mogeneous flow in a cubic box. Unlike the spectra shown in 
Fig. 3 these spectra were computed in the DNS without grav- 
ity. The spectra are vertically shifted for clarity. 


exp-(00/0,)'0"” 


log ;pE(@) a.u. 


-2- 
3 rr ot] 
-| 0 1 2 
l0g;9 M au. 


FIG. 5: The Lagrangian frequency spectra of the particle ve- 
locity fluctuations obtained in a DNS of a particle-laden de- 
caying quantum turbulence. 


Birkhoff-Saffman [19},{13},,[14] 

Ips = ic t)-u(x+r,¢t))dr, (5) 
and Loitsyanskii integrals {15} 

Ip = [Paes t)-u(x+r,¢))dr (6) 
here < ... > denotes the ensemble (or spatial) average. 


Conservation of these integrals by the dissipative 
Navier-Stokes dynamics can be related to the conser- 
vation of the momentum and angular momentum, and 
then to the space homogeneity and isotropy respectively 
(the Noether’s theorem). 


Let us denote a corresponding dynamic invariant as I 
and then let us use this invariant to relate the charac- 
teristic frequency w,. and characteristic velocity u. using 
dimensional considerations 


Ue « Tw (7) 


When the characteristic velocity ue is normally dis- 
tributed the distribution P(w,.) can be readily ob- 
tained from the relationship Eq. (7). Comparison of this 
probability distribution with that provided by the Eq. 
(4) gives a relationship 


2a 


Aes 


(8) 
between the exponents a and £. 


For the Birkhoff-Saffman invariant Eq. (5) we obtain 
6 =1/5and a=3/5. Then from the Eq. (8) one obtains 


E(w) x exp — (w/w) (9) 


For the Loitsyanskii invariant Eq. (6) we obtain 6 = 
1/7 and a = 5/7. Then from the Eq. (8) one obtains 


E(w) x exp —(w/wg)9/2”, (10) 


One can recognize these spectra in the Fig. 1. 


IV. LAGRANGIAN DISTRIBUTED CHAOS IN 
MORE COMPLEX FLOWS 


A. Particle-laden flows 


For the particle-laden flows the randomization can 
be more pronounced in the cases when the gravity is 
not taken into account. Figure 4, for instance, shows 
the Lagrangian frequency spectra obtained in the same 
DNS which was mentioned in Section II (Fig. 3) but 
now without gravity. The spectral data were taken from 
Fig. 13a of the Ref. [10}. One can see that for the corn 
particles 6 < 1 (cf. Fig. 3, the concrete value of 8 = 2/5 
for this case will be discussed in the next sections). For 
the uranium particles the observed value of G = 10/17 
indicates the distributed chaos dominated by the Loit- 
syanskii invariant and corresponding spectrum Eq. (10). 


Figure 5 shows the Lagrangian frequency spectra of 
the particle velocity fluctuations obtained in a DNS 
of a particle-laden decaying quantum turbulence. The 
spectral data were taken from Fig. 10a of the Ref. (16). 


The Gross-Pitaevskii equation was used in this DNS 
and the particles were modeled as potentials having 
classical degrees of freedom described by Newtonian 
dynamics. The well-known Madelung transformation 
was used to transform the quantum description into a 
fluid dynamics-like one. It was already shown (see, for 
instance, Ref. [13}) that the decay of the incompressible 
kinetic energy in the decaying superfluid flow shows 
behavior similar to that observed in decaying classical 
turbulence (mainly because of the sound emission and 
mutual friction effects). 
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FIG. 6: The frequency spectra of the uz component of the 
velocity field measured in the centre of the cubic domain with 
the Rayleigh-Bénard thermal convection. 
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FIG. 7: The Lagrangian frequency spectra of the fluid velocity 
computed in the magnetohydrodynamic DNS of a statistically 
stationary homogeneous isotropic fluid motion for two values 
of the Reynolds number Re = 1050, 3150 (without mean 
magnetic field). The spectra are vertically shifted for clarity. 


The dotted straight line is drawn in the Fig. 5 
for reference to_the Kolmogorov (Lagrangian) scaling 
E(w) «x ew? for the small frequencies (the hard 
- nonsmooth turbulence, see Introduction), and the 
dashed curve indicates the distributed chaos dominated 
by the Loitsyanskii invariant Eq. (10) for large frequen- 
cies (cf Fig. 4). 


B. Rayleigh-Bénard convection 


Let us now apply the above-obtained results to the 
Rayleigh-Bénard thermal convection. The system of 
equations describing the Rayleigh-Bénard convection (in 
the Boussinesq approximation) can be taken as follows 


-Vj)u ers +og0e,+vV7u (11) 
ot po 


| 
+ 
‘e 
oI 
fl 


u-V)O= anil. + «V7, (12) 


V-u=0 (13) 


(see, for instance, Ref. [{19]). Here p and @ are the 
pressure and temperature fluctuation fields respectively 
(9 = T — To(z) here To(z) is a linear vertical profile of 
the temperature), e, is the vertical unit vector, po is the 
mean density of the liquid, g is the gravity acceleration, 
o is the thermal expansion coefficient, v and « are the 
viscosity and thermal diffusivity respectively, AT’ is 
the difference in temperature of the bottom and top 
boundaries, H is the distance between the bottom and 
top boundaries. The periodic boundary conditions are 
usually applied on the side walls, 6 = 0 on the bottom 
and top boundaries as well as the free-slip conditions for 
velocity. 


For the buoyancy-driven dissipative motion the 
Birkhoff-Saffman integral can be generalized in the form 


Ing = [ (u(x,t)-abctr,t)—x 00,4) 6ac+r,t))dr (14) 
as well as the spectrum Eq. (9) (here y = agH/AT). 


Figure 6 shows the frequency spectrum of u, com- 
ponent of the velocity field measured by nine probes 
located at the centre and vertices of a small cube placed 
at the centre of a large cubic domain in a DNS of 
the Rayleigh-Bénard convection (for Rayleigh number 
Ra = 108 and Prandtl number Pr = 1). The spectral 
data were taken from Fig. 4 of the Ref. [20}. Because 
of the absence of a mean velocity at the centre of the 
cubic domain a spatial contamination of the temporal 
spectrum by sweeping effect (see, for instance, bd) and 
references therein) is also absent. Therefore, the mea- 
sured spectrum can be considered as a truly temporal 
one. 


The dashed curve in the Fig. 6 indicates spectrum Eq. 


(9) (the distributed chaos dominated by the generalized 
Birkhoff-Saffman invariant). 


C. Magnetohydrodynamics 


The magnetohydrodynamic equations for incompress- 
ible isotropic homogeneous case in the Alfvénic units are 


1 
OY 2 Tp x KO a8) 
ot Po 


Ob 


— = 2 
Dt V x (ux b) + 7V*b 


(16), 


Vu1=@, vb =0, (17, 18) 
here b = B/,/gop is normalized magnetic field, and f 


is a forcing function. 


It is known that the system (15-18) has the Birkhoff- 
Saffman and Loitsyanskii invariants as well (1.4), (24). 


In the paper Ref. results of a DNS, using the 
Eqs. (15-18) in a three-dimensional spatial box with 
periodic boundary conditions and without a mean 
magnetic field, were reported. In order to maintain a 
statistically stationary system the velocity and magnetic 
field were forced by independent Ornstein-Uhlenbeck 
processes concentrated in a wavenumber shell around 
the wavenumber kr = 3. 


Figure 7 shows the Lagrangian frequency spectra 
of the fluid velocity computed in this DNS for two 
values of the Reynolds number Re = 1050, 3150. The 
spectral data were taken from Fig. 1 of the Ref. 
(three-dimensional case). The dashed curves indicate 
the spectrum Eq. (9) for Re = 3150 and Eq. (10) 
for Re = 1050 (distributed chaos dominated by the 
Birkhoff-Saffman_ and Loitsyanskii invariants respec- 
tively, cf Ref. [23}) . As it was expected the increase of 
the Reynolds number results in a decrease of (3, i.e. the 
randomization becomes stronger with the increase of the 
Reynolds number. 


V. SPONTANEOUS BREAKING OF LOCAL 
REFLECTIONAL SYMMETRY 


The previous consideration of the distributed chaos 
was concentrated on isotropic homogeneous cases with 
global reflectional symmetry. The global reflection sym- 
metry results in zero mean (global) helicity. 

However, spontaneous breaking of the local reflectional 
symmetry can result in the appearance of the nonzero 
point-wise helicity in this case. This phenomenon is an 
intrinsic property of the chaotic/turbulent flows (see, for 
instance, Refs. {25} ,[26}). 

The appearance of the vorticity blobs having nonzero 
helicity and moving with the fluid should accompan, 
this phenomenon (see, for instance, Refs. 7). (31) 
and references therein). These vorticity blobs having 
high relative helicity are a generic property of the 
chaotic/turbulent flows [28). The low viscous dissipation 
over these blobs would result in adiabatic invariance of 


their associated helicity (27), [30]. 


Since the global helicity should be identically equal 
to zero at this phenomenon the localized positive and 
negative helicities associated with the blobs have to be 


canceled at the global average. 


The helicity associated with a vorticity blob is 
Hi; -| h(x, t) dr. (19) 
Vj 


here A(r,t) = u- w is the helicity distribution, V; is the 
spatial volume of the j-blob, and w(r,t) = V x u(r, t) is 
the vorticity field. 


Moments of the helicity distribution h(r,t) = u-w are 
defined as [30] 


; 1 os 
In = Jim 5 oD H! (20) 
j 
where V is the blobs’ total volume. 


Let H; be the total helicity associated with the blobs 
with negative helicity and H fs be the total helicity as- 
sociated with the blobs with positive helicity. Then the 
sign-defined moments of the helicity distribution are 


I+ = lim + SoHE” (21) 


Vow V 


The summation in Eq. (21) is performed for the blobs 
with negative (or positive) Hj only. 


For the odd values of n we have I, = It +I; = 0, due 
to the global reflectional symmetry. Hence, for the odd 
values of n we obtain If = —I;. 


The blobs with high associated helicity should pro- 
vide the main contribution to the moments with high 
values of n (though, the value n = 2 can be often 
considered sufficiently high for the strongly intermittent 
chaotic/turbulent cases). That provides additional 
support for the use of the helicity distribution moments 


I, and [IF as the adiabatic invariants. 


VI. LAGRANGIAN DISTRIBUTED CHAOS 
AND MOMENTS OF HELICITY DISTRIBUTION 


If a flow is dominated by the adiabatic invariant [,, for 
even moments or by the adiabatic invariant [> for odd 
moments we can use the dimensional considerations 


te 0 [En |") wee (21) 


or 


eel ae (22) 

for the even or for the odd moments respectively. For 
the both cases 
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FIG. 8: Lagrangian frequency spectrum of velocity in 
isotropic homogeneous fluid motion: statistically stationary 
(upper curve, Re, = 240) and decaying (lower curve). The 
spectra are vertically shifted for clarity. 
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FIG. 9: Lagrangian frequency spectra of velocity along a pas- 
sive tracer trajectory for Re, = 124 in a laboratory experi- 
ment [34}. The spectra are vertically shifted for clarity. 


Then, from Eq. (8) we obtain 


2(2n — 
ca = _ ~ 
In particular, for the two end cases: 
n >> 1- weakly randomized helical distributed chaos: 
E(k) x exp —(k/kg)'/?, (25) 


and for n = 2 (the Levich-Tsinober invariant [29]) - 
strongly randomized helical distributed chaos: 


E(k) « exp —(k/kg)?/” (26) 
For a moderate randomization n = 3 


E(k) o exp —(k/kg)?/® (27) 


One can recognize spectrum Eq. (27) in the Fig. 1 for 
Rey = 234. This figure shows the Lagrangian frequency 
spectra obtained for statistically stationary isotropic and 
homogeneous chaotic fluid motion computed in a DNS 
[7] (see a description of the DNS in the Introduction). 


An analogous spectrum was obtained in a DNS 
of statistically stationary isotropic and homogeneous 
chaotic motion of an incompressible fluid reported in 
the paper Ref. for Re, = 240. The spectral data 
were taken from Fig. 3 of the Ref. [39}. In this DNS the 
statistically stationary chaotic motion was maintained 
using Lundgren’s linear forcing by applying a force 
proportional to the velocity. 


Figure 8 shows the spectra computed in this DNS 
for two types of isotropic homogeneous motion: statisti- 
cally stationary (upper curve) and freely decaying (lower 
curve). The spectral data were taken from Figs. 1 and 6 
of the Ref. respectively. The dashed curves indicate 
the spectra Eq. (27) for statistically stationary and Eq. 
(10) for decaying cases (cf. Fig. 1 and comments to the 
Fig. 1 in the Introduction). 


VII. LAGRANGIAN HELICAL DISTRIBUTED 
CHAOS IN MORE COMPLEX FLOWS 


One can also recognize spectrum Eq. (27) in Fig. 4 for 
the particle-laden flow (helical distributed chaos for the 
corn particles). If we compare Fig. 4 with Fig. 3 for the 
corn particles we can see that the deterministic chaos 
(8 = 1 for the case with gravitation Fig. 3) has been 
replaced by the helical distributed chaos (3 = 2/5 for the 
case without gravitation Fig. 4), i.e. a randomization 
takes place for the case without gravitation. 


Figure 9 shows the Lagrangian frequency spectra of 
velocity along a passive tracer trajectory for Re, = 124 
in a laboratory experiment (3.4). The spectral data were 
taken from Fig. 5b of the Ref. (3.4). 


In this experiment, the flow of water was generated by 
eight propellers located in the corners of a rectangular 
tank (32 x 32 x 50cm°), and the axis of rotation for 
all propellers was in vertical direction. The rotational 
direction of the propellers was changed in fixed inter- 
vals of time to suppress the existing mean flow. The 
measurements were made by the PTV method using 
neutrally buoyant illuminated particles (with a very low 
Stokes number) inside a spatial ball with a radius of 
50mm centered in the center of the experimental tank. 
Inside this ball the illuminated particles were uniformly 
distributed. 


Figure 9 shows the spectra of the radial and ax- 
isymmetric (vertical) velocity components. The dashed 
curves indicate the spectrum Eq. (27) for both compo- 
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FIG. 10: Frequency spectra of velocity measured with a PIV 
technique for Re, = 220.2; 385 and 555.4 in a laboratory ex- 
periment [35]. 
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FIG. 11: Lagrangian frequency spectra of the fluid veloc- 
ity along the particle trajectory for a heavy-particle-laden 
isotropic homogeneous fluid motion (Re, = 400 and only the 
Stokes drag force was taken into account in the DNS.) 


nents. 


Results of another laboratory experiment with a 
fan-stirred quasi-isotropic homogeneous flow were re- 
ported in paper Ref. (35). The flow was generated 
in a large-scale spherical vessel by four rotating fans 
placed near the wall and was quasi-isotropic and homo- 
geneous in a central spatial ball inside the experimental 
vessel. The velocity characteristics were measured 
with a Particle Image Velocimetry (PIV) technique for 
Rey = 220.2; 385; and 555.4. 


Figure 10 shows the corresponding spectra (the 
spectral data were taken from Fig. 11 of the Ref. ey 
The dashed curves indicate the spectral laws Eq. 27 
for Re, = 220.2 and Eq. (26) for Re, = 385 and 
555.4. As it was expected the value of the dimensionless 
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FIG. 12: Lagrangian frequency spectra of vertical velocity ob- 
tained using the measurements made by a single Lagrangian 
float in the nighttime convective boundary layer. 


parameter @ = 2/5 for Re, = 220.2 is larger then 
that for Re, = 385 and 555.4 (8 = 2/7). Le. the 
randomization of the chaotic flow is lesser for the smaller 
Re, (see Introduction). 


Figure 11 shows the Lagrangian frequency spectra 
of the fluid velocity along the particle trajectory for a 
heavy-particle-laden isotropic homogeneous fluid motion 
(Re, = 400 and only the Stokes drag force was taken 
into account (St is the Stokes number). The DNS 
obtained spectral data were taken from Fig. 6a of the 
Ref. (see also Refs. [37],[38]). The dashed curves 
indicate the spectra Eq. (26) for St = 0.24 and Eq. (27) 
for St = 105. The straight line is drawn for reference 
to the Kolmogorov (Lagrangian) scaling E(w) « ew? 
(the hard - nonsmooth turbulence, see Introduction). 


Interesting examples of the Lagrangian helical dis- 
tributed chaos were obtained in the oceanic subsurface 
observations with the Lagrangian floats [39], [40]. 


In the first example, a subsurface mixed layer (80 m 
deep) was subject to nighttime radiational cooling that 
produced a buoyancy-driven convection. Figure 12 shows 
the Lagrangian spectrum of vertical velocity obtained 
using the measurements made by a single Lagrangian 
float in the nighttime convective boundary layer. The 
spectral data were taken from Fig. 8 of the Ref. [39}. The 
dashed curve indicates the spectrum Eq. (27) (cf Fig. 2). 


In the second example, a Lagrangian frequency 
spectrum was computed from a set of Lagrangian float 
trajectories in Knight Inlet where chaotic/turbulent 
mixing coexists with energetic internal waves. The 
measurements were made by two floats simultaneously 
launched upstream of the Knight Inlet sill in the presence 
of a strong flood tide. 
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FIG. 13: Lagrangian frequency spectra of vertical velocity 


obtained using the measurements made by two Lagrangian 
floats that traverse an inlet sill on strong flood tide. 
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FIG. 14: Lagrangian frequency spectra of the fluid velocity 
obtained in a von Karman-like laboratory ‘swirling’ experi- 
ment (with two counterrotating disks) for Re, = 740. 


Figure 13 shows the Lagrangian frequency spectra of 
the vertical velocity obtained from the float trajectories 
traversing the sill (the_spectral data were taken from 
Fig. 6b of the Ref. [40]). The dashed curve indicates 
the spectrum Eq. (26). 


In a laboratory experiment, described in paper Ref. 
[4i}, a swirling flow was created at high Re, = 740 using 


[1] H. Aref, J. Fluid Mech., 143, 1 (1984) 

[2] U. Frisch and R. Morf, Phys. Rev., 23, 2673 (1981) 

[3] N. Ohtomo, K. Tokiwano, Y. Tanaka et. al., J. Phys. Soc. 
JIpn., 64, 1104 (1995) 

[4] J. E. Maggs and G. J. Morales, Phys. Rev. Lett., 
107,185003 (2011) 

[5] J. E. Maggs and G. J. Morales, Phys. Rev. E 86, 
015401(R) (2012) 


two counterrotating disks (fitted with blades) located at 
the opposite sides of a cylindrical vessel (a von Karman 
type of experimental setting). In the center of the vessel 
the velocity fluctuation can be considered quasi-isotropic 
and homogeneous. Tracking of the tracer particles was 
performed using an acoustic technique. 


Figure 14 shows the power spectrum of the measured 
Lagrangian velocity.The spectral data were taken from 
Fig. 2b of the Ref. [ai]. The dashed curve indicates the 
spectrum Eq. (26). 


VIII. CONCLUSIONS 


The notion of distributed chaos can be used to quan- 
tify the randomization of the Lagrangian chaos. The 
main dimensionless parameter of the distributed chaos - 
8 Eq. (1), has been applied for this purpose. Relevant 
dissipative dynamical invariants - the Birkhoff-Saffman 
and Loitsyanskii integrals, as well as helical dynamical 
invariants (moments of helical distribution) related to 
the spontaneous breaking of local reflectional symmetry, 
have been used to find the corresponding values of ( 
for different flows: isotropic homogeneous fluid motions, 
buoyancy-driven, particle-laden, magnetohydrodynamic, 
chaotic/turbulent motions of quantum fluids, and some 
subsurface oceanic flows (mixing). 


Using this approach it is shown, in particular, that an 
increase of Re, (i.e. intensification of the fluid motion) 
results in an increase of randomization of the Lagrangian 
chaos (see, for instance, Fig. 1, Fig. 7, Fig. 10 and Fig. 
14). For the particle-laden flows gravitation suppresses 
the randomization (cf Fig. 3 with Fig. 4). In some 
cases the Kolmogorov scaling (in the Lagrangian inter- 
pretation E(w) x ew?) can coexist with the Lagrangian 
distributed chaos (in different ranges of scales, see Fig. 5 
and Fig. 11). 
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